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Abstract. Contour dynamics is a classical subject both in physics and in com- 
plex analysis. We show that the dynamics provided by the Lowner-Kufarev ODE 
and PDE possesses a rigid algebraic structure given by the Virasoro algebra. 
Namely, the 'positive' Virasoro generators span the holomorphic part of the com- 
plexified vector bundle over the space of univalent functions, smooth on the bound- 
ary. In the covariant formulation they are conserved by the Lowner-Kufarev evo- 
lution. The 'negative' Virasoro generators span the antiholoniorphic part. They 
contain a conserved term and we give an iterative method to obtain them based on 
the Poisson structure of the Lowner-Kufarev evolution. The Lowner-Kufarev PDE 
provides a distribution of the tangent bundle of non-normalized univalent func- 
tions, which forms the tangent bundle of normalized ones. It also gives an explicit 
correspondence between the latter bundle and the holomorphic eigen space of the 
complexified Lie algebra of vector fields on the unit circle. Finally, we give Hamil- 
tonian and Lagrangian formulations of the motion within the coefficient body in 
the field of an elliptic operator constructed by means of Virasoro generators. We 
also discuss relations between CFT and SLE. 



1. Introduction 

The challenge of structural understanding of non-equilibrium interface dynamics 
has become increasingly important in mathematics and physics. Dynamical inter- 
facial properties, such as fluctuations, nucleation and aggregation, mass and charge 
transport, are often very complex. There exists no single theory or model that can 
predict all such properties. Many physical processes, as well as complex dynamical 
systems, iterations and construction of Lie semigroups with respect to the compo- 
sition operation, lead to the study of growing systems of plane domains. Recently, 
it has become clear that one-parameter expanding evolution families of simply con- 
nected domains in the complex plane in some special models has been governed by 
infinite systems of evolution parameters, conservation laws. This phenomenon re- 
veals a bridge between a non-linear evolution of complex shapes emerged in physical 
problems, dissipative in most of the cases, and exactly solvable models. A sample 
problem is the Laplacian growth, in which the harmonic (Richardson's) moments 
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are conserved under the evolution, see e.g., [231 ES]. The infinite number of evo- 
lution parameters reflects the inflnite number of degrees of freedom of the system, 
and clearly suggests to apply fleld theory methods as a natural tool of study. The 
Virasoro algebra provides a structural background in most of fleld theories, and it is 
not surprising that it appears in soliton-like problems, e.g., KdV or Toda hierarchies, 
see |S1IID|. 

Another group of models, in which the evolution is governed by an inflnite number 
of parameters, can be observed in controllable dynamical systems, where the inflnite 
number of degrees of freedom follows from the inflnite number of driving terms. 
Surprisingly, the same structural background appears again for this group. We 
develop this viewpoint in the present paper. 

One of the general approaches to the growing contour evolution was provided 
by Lowner and Kufarev [201 ES]- The contour evolution is described by a time- 
dependent conformal parametric map from a canonical domain, the unit disk in 
most of the cases, onto the domain bounded by the contour for each flxed instant. 
In fact, these one-parameter conformal maps satisfy the Lowner-Kufarev partial 
differential equation. A characteristic equation to this PDE represents an inflnite 
dimensional controllable system for which the inflnite number of conservation laws 
is given by the Virasoro generators in their covariant form. 

Recently, Friedrich and Werner [8], and independently Bauer and Bernard [1|, 
found relations between SLE (stochastic or Schramm-Lowner evolution) and the 
highest weight representation of the Virasoro algebra. 

All above results encouraged us to conclude that the Virasoro algebra is a common 
structural basis for these and possibly other types of contour dynamics and we 
present the development in this direction here. For the flrst time, a construction, 
which appeared in the fleld theory plays the algebraic structural background for the 
contour evolution in classical complex analysis. 

The structure of the paper is as follows. Sections 2 and 3 contain the necessary 
background on the Virasoro algebra and the Lowner-Kufarev equations. The main 
results are contained in Sections 4 and 5. In Section 4 we construct the Poisson 
structure on the cotangent bundle of the space of univalent functions smooth on the 
boundary and the Hamiltonian system generated by the Lowner-Kufarev equation 
in ordinary derivatives. We establish that the holomorphic Virasoro generators in 
the covariant formulation are conserved under the Lowner-Kufarev evolution (Theo- 
rem 2). The antiholomorphic generators are proved to contain a conserved term and 
we give an iterative method to obtain them based on the Poisson structure of the 
Lowner-Kufarev evolution. The Lowner-Kufarev PDE is shown to provide a distri- 
bution of the tangent bundle of non-normalized univalent functions, which forms the 
tangent bundle of normalized ones. It also gives an explicit correspondence between 
the latter bundle and the holomorphic eigen space of the complexifled Lie algebra 
of vector flelds on the unit circle. In Section 5, we give Hamiltonian and Lagrangian 
formulations of the motion within the coefficient body in the field of an elliptic op- 
erator constructed by means of Virasoro generators. The solutions with constant 
velocity coordinates are found. We prove that the norm of the driving function in the 
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Lowner-Kufarev theory gives the minimal energy of the motion. The short Section 
6 we add for completeness. We briefly review the connections between conformal 
field theory and the Schramm-Lowner evolution following [HIH]. 

Acknowledgements. We are thankful to Helene Airault, Ludwig Faddeev, Paul 
Malliavin, and Yurii Neretin for many helpful discussions concerning the Virasoro 
algebra and its representations. 

2. ViRASORO Algebra 

The Virasoro algebra Vir plays a prominent role in modern mathematical physics, 
both in field theories and solvable models. It appears in physics literature as an 
algebra obeyed by the stress-energy tensor and associated with the conformal group, 
the Virasoro-Bott group, of the worldsheet in two dimensions, see e.g., |2l]. It is a 
unique central extension of the Lie algebra for the Lie-Frechet group Diff of sense- 
preserving diffeomorphisms of the unit circle S"^, and it is an infinite-dimensional real 
vector space. The extension is characterized by a real parameter c, so the Virasoro 
algebra refers to a class of isomorphic Lie algebras corresponding to different values 
of c. At the same time the Virasoro algebra is intrinsically related to the KdV 
canonical structure where the Virasoro brackets become the Magri brackets for the 
Miura transformations of elements of the phase space of the KdV hierarchy (see, 

e.g., mm)- 

The complex hull CVir of the Virasoro algebra can be realized as a central ex- 
tension by C of the Witt algebra, a complex Lie algebra of derivations (or Leibnitz 
rule) of the algebra C[2;,2;~^] of complex Laurent polynomials. The Witt algebra is 
spanned by the generators L„ = -z"^^^ on C\ {0}. The operators L„ plus a central 
element c are called the Virasoro generators. Under any irreducible representation 
of CVir, the quantity c is realized as a complex scalar and is called the central 
charge. The generators satisfy the commutation relations given by 

{Lm,Ln}vir = {n - m)Lm+n + Y^n{n^ - l)Sn-m, {Ln,c}vir = 0, n,meZ, 

where c G C is the central charge. Considering the Virasoro algebra as an operator 
algebra, the generators L„ become the coefficients in a formal Laurent series for the 
analytic component of the stress-energy tensor in 2-D field theory. The attribution 
'Virasoro algebra' is due to a Virasoro's seminal paper [34]. 

Mathematically, the Virasoro algebra appeared for the first time as a central 
extension by the Gelfand-Fuchs cocyde [9] of the Lie algebra Vect of smooth 
vector fields 0^ on the unit circle (see jS]), where the Lie bracket is defined to 
be the commutator of vector fields 

(1) [01,02] = 0102 - 020i. 

Each element of the Lie-Frechet group Diff is represented as 2; = e^"'^^^ with a 
monotone increasing C°° real- valued function a{6), such that a{6 + 2tt) = a{6) + 
2tt. The Lie algebra for this group is identified with Vect S^. The relation of this 
Lie algebra to Diff is subtile because the exponential map is not even locally a 
homeomorphism. 
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2.1. Canonical identification. The entire necessary background of unitary rep- 
resentations of Diff is found in the study of Kirillov's homogeneous Kahlerian 
manifold Diff S^/S^. We deal with the analytic representation of Diff S^/S^. Let 
S stand for the whole class of univalent functions / in the unit disk U normalized 
by f{z) = z{l + Yl^=i (^nZ^) about the origin and C°°-smooth on the boundary 

of U . Given a map / G S we construct the adjoint univalent meromorphic map 

g{z) = diz + dQ + — + . . . , 

z 

defined in the exterior U* = {z : \z\ > 1} of U, and such that C\f{U) = g{U*). Both 
functions are extendable onto 5*^. This conformal welding gives the identification 
of the homogeneous manifold Diff / with the space S: S 9 / ^ o g\gi g 
Diff S'^/S'^, or with the smooth contours F = f{S^) that enclose univalent domains 
fl of conformal radius 1 with respect to the origin and such that oo ^ Q, E Q, see 
[1], [13]. So one can construct complexification of Vect and further projection 
of the holomorphic part to the set A4 C C^, which is the projective limit of the 
coefficient bodies A4 = lim^^oo -A^n, where 

(2) Mn = {{Cl,...,Cn): feS}. 

The holomorphic Virasoro generators can then be realized by the first order differ- 
ential operators 

oo 

Lj = dj + + l)ckdj+k, J G N, 

k=l 

in terms of the affine coordinates of A4, acting over the set of holomorphic functions, 
where dk = d/dck- We explain the details in the next subsection. 

2.2. Complexification. Let us introduce local coordinates on the manifold 
Ai = Diff S^/S^ in the concordance with the local coordinates on the space S 
of univalent functions smooth on the boundary. Observe that Ai is a real infinite- 
dimensional manifold, whereas S is a complex manifold. We are aimed at a complex- 
ification of TAi which admits a holomorphic projection to TS, where Vect qS^ = 
Vect S^/const is a module over the ring of smooth functions, which is associated 
with the tangent bundle TAi. 

Given a real vector space V the complexification Vc is defined as the tensor 
product with the complex numbers V ®r C. Elements of Vc are of the form v ® z. 
In addition, the vector space Vc is a complex vector space that follows by defining 
multiplication by complex numbers, a{v ^ z) = v ^ az for complex a and z and 

V E V. The space V is naturally embedded into V ^ C hj identifying V with 

V <S)l- Conjugation is defined by introducing a canonical conjugation map on Vc as 

V <^ z = V ^ z. 

An almost complex structure J onV can be extended by linearity to the complex 
structure J on Vc by J{v ^ z) = J{v) ® z. Observe that 



J{v ®z) = J{y ® z). 
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Eigenvalues of extended J are ±i, and there are two eigenspaces and V^('''^) 

corresponding to them given by projecting |(1 =F iJ)v. Vc is decomposed into the 
direct sum Vc = V^^'^^ © \/(°'^\ where V^(^'°) = {f O 1 - J{v) ^ i\v G V} and 

are the eigen spaces corresponding to ±i. 

An almost complex structure on Vect O'S'^ may be defined as follows (see [1]). We 
identify Vect qS^ with the functions with vanishing mean value over S^. It gives 

oo 

0(6*) = a„ cos n9 + bn sin n9. 

n=l 

Let us define an almost complex structure by the operator 

oo 

(3) J{(j)){9) = — a„ sin n9 + bn cos n9. 

n=l 

On Vect o^"^ ® C, the operator J diagonalizes and we have the identification 

^ oo 

Vect qS^ 3(l)^v:= -(0 - 2j(0)) = ^(a„ - i6n)e^"^ G (Vect ® C)(^'°), 

^ n=l 

and the latter extends into the unit disk as a holomorphic function. 

The Kirillov infinitesimal action [12] of Vect qS^ on S is given by a variational 
formula due to Schaeffer and Spencer [30l page 32] which lifts the actions from the 
Lie algebra Vect onto S. Let / G S and let 0(e*^) := 4>{9) G Vect be a 
real- valued function in ^ G (0, 27r]. The infinitesimal action 9 ^-^ 9 + e0(e*^) yields a 
variation of the univalent function f*{z) = f + e6vf{z) + o(e), where 

(4) U{z)^^''^'^ f f^f»V <^)d^ 



2n J \ f{w) J w{f{w)-f{z)y 

and ^ f by the above identification. Kirillov and Yuriev [13], [H] (see also [Ij) 
established that the variations 5<^/(C) are closed with respect to the commutator ([T]), 
and the induced Lie algebra is the same as Vect o'S'^. The Schaeffer-Spencer operator 
is linear. 

Treating TAi as a real vector space, the operator 6^ transfers the complex struc- 
ture J from Vect o-S*^ to TA^ by J{S^) := Sj(^). By abuse of notation, we denote 
the new complex structure on TAi by the same character J. Then it splits the 
complexification TA4c into two eigenspaces TAdc = TA4^^'^^ © TAi^^'^\ There- 
fore, 5y = 5^-ij{(j,) ■= 5(j, — iJ{5^) G T^A^^'^\ Observe that 2zdz = —idg on the 
unit circle z = e*^, and = z'^'^^d/dz = —\ie'''^^d/d9 on . Let us take the 
basis of Vect o^*^ (g) C in the form = —ie^^^ in order to keep the index of vec- 
tor fields the same as for L^. Then, the commutator satisfies the Witt relation 
{i^m.,^n} = {n — m)h'n+m- Taking elements Uk = —iw^, \w\ = 1 in the integrand 
of (jlj) we calculate the residue in (jlj) and obtain so called Kirillov operators 

Lj[f]{z)=6.J{z)=z^+'f'{z), J = 1,2,..., 
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SO that these Lj are the holomorphic coordinates on TA4^^''^\ In terms of the affine 
coordinates in A4 we get the Kirillov operators as 

oo 

= + X]*^^ + l)ckdj+k, 

k=l 

where dk = d/dck- They satisfy the Witt commutation relation 

For = we obtain the operator Lq, which corresponds to the constant vectors 
from Vect S^, LQ[f]{z) = zf'{z) — f{z). The elements of the Fourier basis —ie~'^^^ 
with negative indices (corresponding to TA4^^'^^) are extended into U by —iz~''. 
Substituting them in (jlj) we get very complex formulas for which functionally 
depend on Lk (see pLj, [12j), and which are dual to Lk with respect to the action of 
J. The first two operators are calculated as 

L-i[f]iz) = f'{z)~2cj{z)-l, 

L-2[m = ^--^-3ci + (c?-4c2)/(^), 

see [H] . 

This procedure gives a nice links between representations of the Virasoro algebra 
and the theory of univalent functions. The Lowner-Kufarev equations proved to be 
a powerful tool to work with univalent functions (the famous Bieberbach conjec- 
ture was proved [5] using Lowner method). In the following section we show how 
Lowner-Kufarev equations can be used in a representation of the Virasoro algebra. 
In particular, we identify TAi^^''^^ with TAi, equipped with its natural complex 
structure given by coefficients of univalent functions, by means the Lowner-Kufarev 
PDF. 

3. Lowner-Kufarev Equations 

A time-parameter family Q{t) of simply connected hyperbolic univalent domains 
forms a Lowner subordination chain in the complex plane C, for < t < r (where 
T may be oo), if Q{t) ^ ^{s), whenever t < s. We suppose that the origin is an 
interior point of the Caratheodory kernel of {0(t)}[^o- 

A Lowner subordination chain Q{t) is described by a time-dependent family of 
conformal maps z = f{C,t) from the unit disk U = {( : \C\ < 1} onto ^l{t), 
normalized by /(C,^) = oi(t)C + a2(t)C^ + . . . , ai(t) > 0, di(t) > 0. After Lowner's 
1923 seminal paper [20] a fundamental contribution to the theory of Lowner chains 
was made by Pommerenke [251 E6j who described governing evolution equations in 
partial and ordinary derivatives, known now as the Lowner-Kufarev equations due 
to Kufarev's work [TB] . 

One can normalize the growth of evolution of a subordination chain by the con- 
formal radius of Q{t) with respect to the origin by ai{t) = e*. 

Lowner [20] studied a time-parameter semigroup of conformal one-slit maps of 
the unit disk U arriving then at an evolution equation called after him. His main 



VIRASORO ALGEBRA IN LOWNER-KUFAREV DYNAMICS 



7 



achievement was an infinitesimal description of tlie semi-flow of such maps by the 
Schwarz kernel that led him to the Lowner equation. This crucial result was then 
generalized in several ways (see [26] and the references therein). 

We say that the function p is from the Caratheodory class if it is analytic in U, 
normalized as p{() = 1 + pi( + P2C^ + • • • , ( ^ U, and such that Re p{() > 
in U. Pommerenke [25l[26] proved that given a subordination chain of domains Q{t) 
defined for t G [0, r), there exists a function p{(, t), measurable in t G [0, r) for any 
fixed z E U, and from the Caratheodory class for almost all t G [0, r), such that the 
conformal mapping f : U —>■ Q{t) solves the equation 

for C & U and for almost all t G [0, r). The equation ([5]) is called the Lowner-Kufarev 
equation due to two seminal papers: by Lowner [20] who considered the case when 

iu{t) I /- 

(6) P(C^^) = ^m^^^ 

where u{t) is a continuous function regarding to t G [0, r), and by Kufarev [16] who 
proved differentiability of / in t for all ( from the kernel of {fl{t)} in the case of 
general p in the Caratheodory class. 

Let us consider a reverse process. We are given an initial domain f2(0) = Qq (and 
therefore, the initial mapping /(C^O) = /o(C))! and a function p(C,^) of positive 
real part normalized by ^(C,^) = 1 + PiC + • • • • Let us solve the equation (jS]) and 
ask ourselves, whether the solution f{C,t) defines a subordination chain of simply 
connected univalent domains f{U,t). The initial condition /(C^O) = /o(C) is not 
given on the characteristics of the partial differential equation ([5]), hence the solution 
exists and is unique but not necessarily univalent. Assuming s as a parameter along 
the characteristics we have 

dt dC , . df 

as as as 

with the initial conditions t(0) = 0, ({0) = z, f{(,0) = fo{(), where z is in U. 
Obviously, t = s. Observe that the domain of ( is the entire unit disk. However, the 
solutions to the second equation of the characteristic system range within the unit 
disk but do not fill it. Therefore, introducing another letter w (in order to distinguish 
the function w{z,t) from the variable () we arrive at the Cauchy problem for the 
Lowner-Kufarev equation in ordinary derivatives 

/^\ dw . , 

(7) — = -wp{w,t), 

for a function ( = w{z, t) with the initial condition w{z, 0) = z. The equation (jTj) is a 
non-trivial characteristic equation for (JSj). Unfortunately, this approach requires the 
extension of /o(ui~^(C, t)) into the whole U {w~^ means the inverse function) because 
the solution to ([5]) is the function f{C,t) given as fo{w~^{(,t)), where ( = w{z, s) is 
a solution of the initial value problem for the characteristic equation ([7]) that maps 
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U into U. Therefore, the solution of the initial value problem for the equation ([5]) 
may be non- univalent. 

Solutions to the equation ([7j) are holomorphic univalent functions w{z, t) = e~^z + 
a2{t)z'^ + . . . in the unit disk that map U into itself. Every function / from the class 
S can be represented by the limit 



where w{z,t) is a solution to ([7]) with some function p{z,t) of positive real part 
for almost all t > (see [26], pages 159-163]). Each function p{z,t) generates a 
unique function from the class S. The reciprocal statement is not true. In general, 
a function / G S can be obtained using different functions p{-,t). 

Now we are ready to formulate the condition of univalence of the solution to the 
equation ([5]), which can be obtained by combination of known results of [26] . 

Theorem 1. [26l [27| Given a function p{C,t) of positive real part normalized by 
P(C; t) = 1 + pi( + . . . , the solution to the equation ^ is unique, analytic and 
univalent with respect to ( for almost all t > 0, if and only if, the initial condition 
/o(C) is taken in the form (jS]), where the function w{(,t) is the solution to the 
equation ([7]) with the same driving function p. 

Recently, we started to look at Lowner-Kufarev equations from the point of view 
of motion in the space of univalent functions where Hamiltonian and Lagrangian 
formalisms play a central role (see, [32]). Some connections with the Virasoro al- 
gebra were also observed in [22l[32]. The present paper generalizes these attempts 
and gives their closed form. The main conclusion is that the Lowner-Kufarev equa- 
tions are naturally linked to the holomorphic part of the Virasoro algebra. Taking 
holomorphic Virasoro generators L„ as a basis of the tangent space to the coefficient 
body for univalent functions at a fixed point, we see that the driving function in 
the Lowner-Kufarev theory generates generalized moments for motions within the 
space of univalent functions. Its norm represents the energy of this motion. The 
holomorphic Virasoro generators in their co-tangent form will become conserved 
quantities of the Lowner-Kufarev ODE. The Lowner-Kufarev PDE becomes a tran- 
sition formula from the affine basis to Kirillov's basis of the holomorphic part of 
the complexified tangent space to Ai at any point. Finally, we propose to study an 
alternate Lowner-Kufarev evolution instead of subordination. 

4. Witt algebra and the classical Lowner-Kufarev equations 

In the following subsections we reveal the structural role of the Witt algebra as 
a background of the classical Lowner-Kufarev contour evolution. As we see further, 
the conformal anomaly and the Virasoro algebra appear as a quantum or stochastic 
effect in SLE. 

4.1. Lowner-Kufarev ODE. Let us consider the functions 



f{z) = lim e^w{z, t), 
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satisfying the Lowner-Kufarev ODE 

dw , , 

(9) — = -wp{w,t), 

with the initial condition w{z, 0) = z, and with the function p{z, t) = 1 +pi{t)z + . . . 
which is holomorphic in U and measurable with respect to t G [0, oo), such that 
Re p > in f/. The function w{z,t) is univalent and maps U into U. 

Lemma 1. Let the function w{z,t) be a solution to the Cauchy problem for the 
equation ([9]) with the initial condition w{z,0) = z. If the driving function p{-,t) , 
being from the Caratheodory class for almost all t > 0, is C°° smooth in the closure 
U of the unit disk U and summable with respect to t, then the boundaries of the 
domains B(t) = w{U, t) dU are smooth for all t. 

Proof. Observe that the continuous and differentiable dependence of the solution to 
a differential equation x = F{t, x) on the initial condition x{0) — Xq is cL classical 
problem. One can refer, e.g., to [55] in order to assure that summability of F{-,x) 
regarding to t for each fixed x and continuous differentiability {C^ with respect to 
X for almost all t) imply that the solution x{t, Xq) exists, is unique, and is with 
respect to xq. In our case, the solution to IQ exists, is unique, analytic in U, and 
moreover, on its boundary S^. Let us differentiate iQ inside the unit disk U 
with respect to z and write 

t 

logw' = - {p{w{z,t),t) + w{z,T)p'{w{z,r),T))dT, 



choosing the branch of the logarithm such as logw'{0,t) = —t. This equality is 
extendable onto because the right-hand side is, and therefore, w' is and w is 
on S^. We continue analogously and write the formula 

t 

w" = —w' J {2w'{z, t)p'(w{z, t),t) + w{z, t)w'{z, t)p"{w{z, r), T))dT, 



which guarantees that w is on S^. Finally, we come to the conclusion that w is 
C°° on S\ □ 

Let f{z,t) denote e^w{z,t). The limit lim^^oo /(-z, ^) is known [25] to be a repre- 
sentation of all univalent functions. 

Let the driving term p{z., t) in the Lowner-Kufarev ODE be from the Caratheodory 
class for almost alH > 0, C°° smooth in f/, and summable with respect to t. Then 
the domains VL{t) = w{U,t) have smooth boundary dQ{t). So the Lowner equation 
can be extended onto the closed unit disk U = U U S^. 

Consider the Hamiltonian function given by 

r - dz 

(10) H= / f{z,t){l-p{e-'f{z,t),tmz,t)-, 

J iz 
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on the unit circle z G S"^, where ip{z,t) is a formal series 

oo 
n=—k 

defined about the unit circle 5*^ for any A; > 0. The Poisson structure on the 
symplectic space (/, tp) is given by the canonical brackets 

or in coordinate form (only ijjn for n > 1 are independent co- vectors corresponding to 
the tangent vectors dn with respect to the canonical Hermitean product for analytic 
functions) 

, , dp dq dp dq 

Here 

P{t)= [ p{z,t)-, Q{t)= [ q{z,t)-. 

J IZ J IZ 

The Hamiltonian system becomes 

(11) = /(I - P(e-V, ^)) = ^ = {/, H}, 

at Sip 

for the position coordinates and 

(12) ^ = - t) - e-'fp'{e-'f, m = ^ 

for the momenta, where ^ and 4f are the variational derivatives. So the phase 

_ (>i Sip ^ 

coordinates (/, ip) play the role of the canonical Hamiltonian pair. 

The coefficients c„ are the complex local coordinates on A^, so in these coordinates 
we have 

dz 



dc 6^ f 

= Cn - — / W{z,t)p{w{z,t),t) 

1 rJ^ dz 



I Y^e-'-(e'wf»p,^,, „>1. 

i, k=l 



2Txi 

s 

Let us fix some n and project the infinite dimensional Hamiltonian system on an 
ra-dimensional A^„. The dynamical equations for momenta governed by the Hamil- 
tonian function (fTOl) are 



1 ^ f dz 

= -^3 + TT" / + ^pO k j+i ' i = i> • • • ) ^ - 1) 

1 — 1 ^ 



and 

(13) = 0. 
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In particular, 

Cl 
C2 
C3 

for ri = 3 we have 

^1 

Let us set the function L{z) := f'{z,t)'ijj{z,t). Let {L{z))^o mean the part of the 
Laurent series for L{z) with negative powers of z, 

1 - - 1 °° L 

(L(2;))<o = (^1 + 2ci^2 + 3c2V^3 + ...)- + (V^2 + 2ciV^3 + • • • )^ + ■ ■ ■ = 5^ ^• 

fc=l 

Then, the functions i^(-z) and {L{z))^q are time-independent for all z G S^. 

It is easily seen that, passing from the cotangent vectors ipk to the tangent vectors 
dk, the coefficients Lk of {L{z))^o defined on the tangent bundle TAi^^'^^ are exactly 
the Kirillov vector fields Lk- The corresponding fields Lk in the covariant form are 
conserved by the Lowner-Kufarev ODE because = {L^, H} = 0. The above 
Poisson structure coincides with that given by the Witt brackets introduced for 
previously. For finite-dimensional grades this result was obtained in |22j . 

Let us formulate the result theorem. 

Theorem 2. Let the driving term p{z,t) in the Lowner-Kufarev ODE be from the 
Caratheodory class for almost all t > 0, smooth in U , and summahle with respect 
to t. Then the Kirillov fields in the covariant form are the conserved quantities for 
the Hamiltonian system [IJIW^) generated by the Lowner-Kufarev ODE. 

Remark 1. Another way to construct a Hamiltonian system could be based on the 
symplectic structure given by the Kdhlerian form on Diff S^/S^. However, there is 
no explicit expression for such form in terms of functions / G S. Moreover, there 
must be a Hamiltonian formulation in which the Lowner-Kufarev equation becomes 
an evolution equation. This remains an open problem. 

Remark 2. At a first glance the situation with an ODE with a parameter is quite 
simple. Indeed, if we solve an equation of type f{t,e^^) = F{f{t,e^^),t), then fixing 6 
we have an integral of motion C = I{f{t, ■),{) = const. Then, releasing 9, we have 
C(e*^) = I{f{t,e^^),t). Expanding C(e*^) into the Fourier series, we obtain an 
infinite number of conserved quantities, but they do not manifest an infinite number 
of degrees of freedom that govern the motion as in the field theory where the governing 
equations are PDE. In our case, we have not only one trajectory fixing the initial 
condition but a pensil of trajectories because our equation has an infinite number 
of control parameters, the Taylor coefficients of the function p{z,t), which form a 



= -e~'pi, 

= -2e"ViCi - 6-2^2, 

= -e-'pi{2c2 + cl) - 3e-^'p2Ci - e'^Vs, 

= 2e-Vi^2 + (2e-ViCi + 3e->2)V^3, 
= 2e-Vi^3, 
= 0. 
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bounded non-linear set of admissible controls. Therefore, we operate with sections of 
the tangent and co-tangent bundles to the infinite dimensional manifold A4. instead 
of vector fields along one trajectory as in usual ODE. 

Remark 3. No linear combinations of Li, . . . , Ln, . . . allows us to reduce the 
system of {Lk} to a new system of involutory {^^} in order to claim the Liouville 
integrability of our system. Observe that the coefficients in these linear combinations 
must be constants to keep conservation laws. 

4.2. Construction of Lq and L_„. Consider again the generating function L{z) — 
f'{z,t)'4){z,t) and the 'non-negative' part {L{z))>q of the Laurent series for L{z), 

{L{z))>o = (V'O + 2CiVSl + 302^^2 + ...) + (V'-l + SCiVSq + 302^^1 + ...)z + ... 

oo 

k=0 

All are conserved by the construction. Define ipo = — Yl'^=i Ck'ipk, and 

Lo = £o - (^0 - ^o)- 

The operator Lq acts on the class S by LQ[f]{z) = zf'{z) — f{z). Next define 
L_i = C_i - - 2p*_^) - 2ci{i!o - where V-* i = 0. Then, 

L_^[f]{z) = f'{z)-2c^f{z)-l 

Finally, 

L-2 = >C_2 - (V'-2 - - 2Cl('0-l - - 3C2('0o " 

We choose ■0I2 — (c3 — 3ciC2 + cf + . . . , so that 



and 



r-2[m = \-^^-c,-{c2-c^)f{z), 

LArnz) = - - 3ci + (c? - 4c2)/(z). 



An important fact is that 

Lq = Citpi + 2C2V^2 + • . ■ , 

L„i = (3c2 - 2cl)t/ji + . . . , 
L-2 = (5c3 - 6cic2 + 2cl)^i + ..., 

are linear with respect to ipk, ^ > 1, and therefore, are sections of T*A4, which are 
dual to Kirillov's vector fields. Equivalently, 

Lo -1 _2[/](^) = function(ci, C2, . . . + . . . , 2;*^ = . 

All other co-vectors we construct by our Poisson brackets as 

1 1 _ 1 r 
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The form of the Poisson brackets guarantees us that all L_„ are linear with respect 
to %jji,tjj2, . . . and span the anti-holomorphic part of the co-tangent bundle T^'^'^^ Ai. 

Let us summarize the above in the following conclusion. We considered a non- 
linear contour dynamics given by the Lowner-Kufarev equation. It turned out to be 
underlined by an algebraic structure, namely, by the Witt algebra spanned by the 
Virasoro generators L„, G Z. 

• L„, n = 1,2, .. . are the holomorphic Virasoro generators. They span the 
holomorphic part of the complexied tangent bundle over the space of univa- 
lent functions, smooth on the boundary. In the covariant formulation they 
are conserved by the Lowner-Kufarev evolution. 

• Lq is the central element. 

• L^n, n = 1,2, . . . are the antiholomorphic Virasoro generators. They span 
the antiholomorphic part of the decomposition. They contain a conserved 
term and we give an iterative method to obtain them based on the Poisson 
structure of the Lowner-Kufarev evolution. 

4.3. Lowner-Kufarev PDE. The Lowner equation in partial derivatives is 
^(C, t) = Cw'iC, t)p{C, t). Re p(C, t) > 0, Id < 1. 

with some initial condition w{z,0) = fo{z). Let us consider the one-parameter 
family of functions f{z,t) = e^^w{z,t) = z{l + J2'^=i^n(t)z"'), f{z,0) = fo{z) as a 
path in S. At the initial point foiz) we have that Tf^S = Tf^M^^'^^ = Tf^M. A 
path in the coefficient body A4 in the neighbourhood of /o is (ci(t), . . . , c„(t), . . . ) 
with the velocity vector Cidi + ■ ■ ■ + c„9„ + ■ ■ ■ G Tf^Ai. 

Taking the Virasoro generators {L^}, A; > 1, as a basis in T^Ai^^'^^ we wish the 
velocity vector written in this new basis to be 

(14) cidi H h Cndn H = UiLi + ... UnL^ + .... 

We compare f|T^ with the Lowner-Kufarev equation 

(15) / = Cidi H h Cndn H = Zf'p{z, t) - / = Lo + UiLi + ... UnL^ + 

where p{z, t) = 1 + uiz + ■ ■ ■ + Unz" + . . . , and Lq/ = zf — /. In view of similarity 
between these two expressions f fT^ and f|T5|) . we notice that 

• a new term Lq appears in the Lowner-Kufarev equation; 

• the function p{z, t) with positive real part corresponds to subordination, 
whereas for generic trajectories it may have real part of arbitrary sign. We 
call this an alternate Lowner-Kufarev evolution; 

• the vector Lq corresponds exactly to the rotation: 

e-/(e-z) = f{z) - ze{zf'{z) - f{z)) + o{e). 

Let us consider the set Sq of non-normalized smooth univalent functions of the 
form F{z, t) = ao{t)z + ai{t)z'^ + . . . , with a tangent vector doc^o + ■ — I- o,ndn + • • • , 
where dk = d/dak, k = 0,1,2, ... . Our aim is to define two different distributions 
for the tangent bundle TSq, that form a sub-bundle of co-dimension 1, which is the 
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tangent bundle TS. This will be realized by means of formulas (fT^ and (fT5|) . Notice 
that dkF = z^+^. Setting Lfc(F) := z'^+'^F' we get 

F = dodo H h dndn H = zf'p{z, t) = uoLq + uiLi + . . . m„L„ + . . . , 

where t) = uq + Mi^; + ■ ■ ■ + UnZ^ + . . . . This alternate Lowner-Kufarev equation 
represents recalculation of the tangent vector in the new basis 

oo^o H \- «n<9„ H = uoLo + . . . m„L„ + . . ., 

where Lk = aodk + 2aidk+i + 

Let us present the distributions. We start with F G Sq, then we define / G S. 
The necessary distribution is the map 

So3F^TfS^ TpSo. 

The analytic form of the first distribution is the following factorization fi{z,t) = 
-F(z, t) = z + + . . so that 

ao ^ ' ^ ao ' 

(16) h = zf[p{z,t)-^h, 

ao 

where Uo = —. Then we obtain 

cidi H h Cndn H = Lo + uiLi H h m„I/„ + • • • 

where Lq/i = uo{zf[ - fi), Lkfi = z''+^f[, Ck = ^, dk = In particular, oq = e* 
implies the Lowner-Kufarev equation for arbitrary sign of Re p. 

The analytic form of the second distribution becomes f2iz,t) = F{-^z,t) = z + 
^z^ + .... so that 

(17) h = zfyi-,t)-^zf^, 

ao ao 

where again mq = ^. In the coefficient form we get 

Cidi H h Cndn H = UiLi H h m„L„ + • • • 

where Lfc/a = z^-^^f^, ct = dk = 

Observe that the equation (fT7|) gives an identification of TAi^^'^^ with TA^. 
Finally, let us make an explicit calculation of Lo, which for ao = e* we continue 
to denote by Lq. Using Kirilov's basis Li, L2, ... as a linear combination we write 

00 

Lo = ^ ] ^rriLm- 

m=l 

The coefficients 11^ are polynomials, which can be obtained using the following 
recurrent formulas 

m— 1 m 
Kl = 0, K^ = - ^j{m- j + l)Cm-jCj, Um = mCm + ^ Km-j+lPj-1, 
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where Pk are polynomials 

k 

(18) Po = l, Pi = -2ci, P2=Acl-3c2, Pfc = -5^(j + l)c,Pfc_„ 

i=i 

Let us summarize the above considerations in the following theorem. 

Theorem 3. The Lowner-Kufarev PDF flTBl) gives the distribution for the tangent 
bundle TSq of non-normalized smooth univalent functions So, that forms a sub- 
bundle of CO- dimension 1, which is the tangent bundle TS. 

The equation (fT7|) gives another distribution, and moreover, it makes the explicit 
correspondence between the natural complex structure ofTS, as S embedded into , 
and the complex structure ofTAi^^''^^ at each point f E S defined by ([3]). 

One of the reason to consider the alternate Lowner-Kufarev PDE is the regularized 
canonical Brownian motion on smooth Jordan curves. For all Sobolev metrics if 2+<^, 
the classical theory of stochastic flows allows to construct Brownian motions on 
diffeomorphism group of S^. The case 3/ 2 is critical. Malliavin [21] constructed the 
canonical Brownian motion on the Lie algebra Vect for the Sobolev norm 
Another construction was proposed in |7j. Airault and Ren [2] proved that the 
infinitesimal version of the Brownian flow is Holder continuous with any exponent 
P <1. 

The regularized canonical Brownian motion on Diff 5*^ is a stochastic flow on 
associated to the Ito stochastic differential equation 

, „ {,X2n{t) cosnO - X2n~i{t) sinra6'), 

n=l * 

where {xk} is a sequence of independent real-valued Brownian motions and r e (0, 1) 
and the series for (xti^) ^ Gaussian trigonometric series. Kunita's theory of 
stochastic flows asserts that the mapping 6 —>■ gl.i{9) is a diffeomorphism and 
the limit lim g^^ j = g^^t exists uniformly in 9. The random homeomorphism g^^t 

is called canonical Brownian motion on Diff S^, see [21 [71 [211 12H]- If was shown in 
[21 [7|, that this random homeomorphism is Holder continuous. 

The canonical Brownian motion can be defined not only on Diff S^, but also on the 
space of C°°-smooth Jordan curves by conformal welding. This leads to dynamics 
of random loops which are not subordinated. 

5. Elliptic operators over the coefficient body 

The Kirillov first order differential operators generate the elliptic operator 
|i^fcp. In this section we construct the geodesic equation and find geodesies with 
constant velocity coordinates in the field of this operator. In particular, we shall 
prove that the norm of the driving function in the Lowner-Kufarev theory gives the 
minimal energy of the motion in this field. 
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5.1. Dynamics within the coefficient body. Let us recall the geometry of the 
coefficient body A4n for finite n. The affine coordinates are introduced by projecting 

oo 

M3 f = z{l + ^CkZ^'^ ^ (Ci,. . . ,C„) G Mn. 
k=l 

The manifold Ain was studied actively in the middle of the last century, see e.g., 
[31 [30]. We compile some important properties of A1„ below: 

(i) Ain is homeomorphic to a (2n — 2)-dimensional ball and its boundary dAin 
is homeomorphic to a (2n — 3)-dimensional sphere; 

(ii) every point x G dAin corresponds to exactly one function / G S which is 
called a boundary function for Ain', 

(iii) boundary functions map the unit disk U onto the complex plane C minus 
piecewise analytic Jordan arcs forming a tree with a root at infinity and 
having at most n tips, 

(iv) with the exception for a set of smaller dimension, at every point x G dAin 
there exists a normal vector satisfying the Lipschitz condition; 

(v) there exists a connected open set Xi on dAin, such that the boundary dAin 
is an analytic hypersurface at every point of Xi. The points of dAin cor- 
responding to the functions that give the extremum to a linear functional 
belong to the closure of Xi. 

Properties (ii) and (iii) imply that the functions from S deliver interior points of 
A4n- The Kirillov operators Lj restricted onto Atn give truncated vector fields 

n-j 

= + + l)ckdj+k, 

k=l 

which we, if it causes no confusion, continue denoting by Lj in this section. In [22] 
based on the Lowner-Kufarev representation, we showed that these Lj can be ob- 
tained from a partially integrable Hamiltonian system for the coefficients in which 
the first integrals coincide with Lj. 

Let c{t) = (ci(t), . . . , Cn(t)) be a smooth trajectory in Ain', that is a map 
c : [0, 1] Ain- Then the velocity vector c(t) written in the affine basis as c{t) = 
Ci(t) di + . . . + Cn(t) dn can be also represented in the basis of vector fields Li, . . . , L„ 
(compare with (ITT]) ) as 

(19) C{t) =Ci{t)di + ... + Cn{t)dn 

= UiLi + U2L3 + . . . + UnLn, 

where the coefficients Uk can be written in the recurrent form as 

fc-i 

(20) ui = ci, Uk = Ck- ^(j + l)cjUk-j. 
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Expressing Uk in terms of and c^, we get 

fc-i 

(21) Uk = Ck + PjCk-j- 

i=i 

One may notice that these polynomials are the first coefficients of the holomorphic 
function l/f'{z), where / G S. In the infinite dimensional case this follows from 
the Lowner-Kufarev equation (fTTl) with oq = e*. Kirillov's fields Lk act over these 
polynomials as 

LkPn = {n — 2k — l)Pn-k n > k and L^Pn = n < k. 
Proposition 1. We define 

uji = dci, 

002 = dC2 — 2CiUJi, 



n-1 

(22) UJn = dCn - X^(j + l)CjUJn-j. 

i=i 

Then, {ui, . . . , Un} is a conjugate to {Li, . . . , L„} basis of one-forms. Namely, 

LUniLn) = 1, t^n(-^fc) =0 if k 71. 

Proof. If > n, then the vector fields do not contain 9„. Since the form uOn 
depends only on dcj with j < n, then 

n-1 

uJn{Lk) = dn{Lk) - ^(j + l)cjUJn-j{Lk) =0 for k > n > n - J. 

3=1 

If n = k, then 

n-1 

UJn{Ln) = dn{Ln) - ^(j + l)CjUn-j{Ln) = 1 + for n> H-j. 

i=i 

To prove the case k < n we apply the induction. Let us show for Li. We have 

uj2iLi) = dc2{Li) - 2ci(Li) = 2ci - 2ci = 0. 
We suppose that ci;„(Li) = 0. Then 

n 

UJn+l{Li) = dCn+l{Li) - ^(j + l)CjUJn+l-j{Li) = {u + l)c„ - {u + l)c„t^i(Li) = 0. 

i=i 

The same arguments work for uJn{Lk) with k < n. □ 

In the affine basis the forms uJk can be written making use of the polynomials 
Pn- We observe that one-forms are defined in a similar way as the coordinates 
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Uk with respect to the Kirillov vector fields L^. Thus, if we develop the recurrent 
relations (!22|) and collect the terms with dcn we get 

k-l 

Uk = dck + ^ Pjdck-j. k = 1, . . . ,n. 
i=i 

By the duality of tangent and co-tangent bundles the information about the mo- 
tion is encoded by these one-forms. 

5.2. Hamiltonian equations. There exists an Hermitian form on TAi^, such that 
the system {Li, . . . , L„} is orthonormal with respect to this form. The operator L = 
l-^fcP is elliptic, and we write the Hamiltonian function H {c, c, ip , ip) defined on 
the co-tangent bundle, corresponding to the operator L as H{c, c, ip, ip) = "^l^i 
where 

n—k 

ik = 'ipk + ^(j + i)cji!k+j- 

The corresponding Hamiltonian system admits the form 

dH - 

Ci = = li 



In + ^(j + l)Cjln-j 



n—p 



ip = = + i)X^4^fc+p 



P k = l 



Vn = -T. = 0. 

OCn 



Let us observe that 



n—k 



(23) ik = ^(j - k)ljlj+k- 
Expressing Ik from the first n Hamiltonian equations we get 

k-l 

(24) 4 = Cfc + ^ PjCk-j, k = l,...,n. 

i=i 

We can decouple the Hamiltonian system making use of (!23|) and (12^ which leads 
us to the following non-linear differential equations of the second order 

k-l k-l 

Ck = h + 

i=i i=i 
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where Ik are expressed in terms of the product of Ijlj+k by fl23|) . and the last products 
depend on Pj, Pj and c, Cj for the corresponding indices j by ( l24l) . For example, 



j=i p=i 



pCj^p 



5=1 



Comparing and (12T!) . we conclude that Ik = Uk and Uk satisfy the differential 
equations 



(25) 



n—k 

ilk = ^(j - k)ujUj+k, 



on the solution of the Hamiltonian system. Observe that any solution of ( l25i) has a 
velocity vector of constant length. It is easy to see from the following system 



UiUi 

U2U2 

(26) 

U4U4 

Then, 



OU1U1U2 + U1U2US + 2uiu^U4 + 3U1U4U5 + Auiu^ue + . . . 

— IU1U2U3 + OU2U2U4 + '^U2U3U5 + 2U2U4UQ + . . . , 

— 2UiUsU4 — IU2U3U5 + OUsU^Uq + . . . , 

—3UIU4U5 — 2U2U4UQ + . . . , 

-AuiU^Uq + . . . , 



d\u\- 



2 '^iukUk + UkUk) 

k=l 



0, 



for any n, thanks to the cut form of our vector fields and the skew symmetry of (1261) . 
The simplest solution may be deduced for constant driving terms Uk, k = 1, . . . ,n. 
The Hamiltonian system immediately gives the geodesic 



Ci = Mi(0)s + Ci(0), 

C2 = M?(0)s2 + M2(0)s + C2(0), 



C3 



3 2 2 

3mi(0) {uliOYj + U2{0yj + C2(0)) + 2M2(0)(Mi(0)y + Ci(O)s) + Us{0)s + 02(0), 



In general, c„ becomes a polynomial of order n with coefficients that depend on the 
initial data c(0) and on the initial velocities u{0). 

The Lagrangian C corresponding to the Hamiltonian function H can be defined 
by the Legendre transform as 

n k—1 ^ n 

C = {c,ij)-H = J2 {h^k + ^k Y^{3 + l)c,-4-,) - 2 5Z l^'^l'- 

k=l j=l k=l 
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Taking into account that 



fc-i 



Summing up over /c, we obtain (c, ip) = X]fc=i h^k = Sfc=i ""fc^fc, that gives us 



2 

fc=i 



All these considerations can be generalized for n oo. Thus, we conclude that 
the coefficients of the function p{z, t) in the Lowner-Kufarev PDE play the role of 
generalized moments for the dynamics in and M. with respect to the Kirillov 
basis on the tangent bundle. Moreover, the L^-norm of the function p on the circle 
is the energy of such motion. 



6. SLE AND CFT 

In this section we briefly review for completeness the connections between con- 
formal field theory (CFT) and Schramm-Lowner evolution (SLE) following, e.g., 
[HIH]). SLE (being, e.g., a continuous limit of CFT's archetypical Ising model at its 
critical point) gives an approach to CFT which emphasizes CFT's roots in statistical 
physics. 

SLE^ is a x-parameter family of covariant processes describing the evolution of 
random sets called the SLE^ hulls. For different values of x these sets can be either 
a simple fractal curve x G [0, 4], or a self-touching curve x G (4, 8), or a space ffiling 
Peano curve x > 8. At this step we deal with the chordal version of SLE. The 
complement to a SLE^ hull in the upper half-plane EI is a simply connected domain 
that is mapped conformally onto H by a holomorphic function g{z,t) satisfying the 
equation 

dt g{z,t)-^t 

where = y/^Bt, and Bt is a normalized Brownian motion with the diffusion 

2t 

constant x. The function g{z, t) is expanded as g{z,t) = zA V ■ ■ ■ ■ The equation 

(j27|) is called the Schramm-Lowner equation and was studied first in |T7]-|19], see 
also [29] for basic properties of SLE. Special values of x correspond to interesting 
special cases of SLE, for example x = 2 corresponds to the loop-erasing random 
walk and the uniform spanning tree, x = 4 corresponds to the harmonic explorer 
and the Gaussian free field. Observe, that the equation ( 1271) is not a stochastic 
differential equation (SDE). To rewrite it in a stochastic way (following |3], [S]) let 
us set a function kt(z) = g{z,t) — C,t, where kt{z) satisfies already the SDE 

2 

dkt{z) = -r-r-rdt - d^f 
kt{z) 
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For a function F{z) defined in the upper half-plane one can derive the Ito differential 
(28) dF{kt) = -d^tL^iFih) + dt{^L\ - 2L^2)F{kt), 

with the operators L_i = — ^ and L_2 = —jj^- These operators are the first two 
Virasoro generators in the 'negative' part of the Witt algebra spanned by the oper- 
ators —z"''^^-^ acting on the appropriate representation space. All other generators 
can be obtained by the commutation relation 

For any state {ip), the state L_i\ip) measures the diffusion of \ip) under SLE, and 
{^L'^i — 2L_2)\'ip) measures the drift. The states of interest are drift-less, i.e., the 
second term in (1281) vanishes. Such states are annihilated by ^L'^i — 2L_2, which 
is true if we choose the state as the highest weight vector in the highest weight 
representation of the Virasoro algebra with the central charge c and the conformal 
weight h given by 

(6-x)(3k-8) , 6-K 

c = , h = , 

2x 2x 

and the operators L_i and L_2 are taken in the corresponding representation. It was 
obtained in [1] and [8], that F{kt) is a martingale if and only if {^L'^^ — 2L_2)F{kt) = 
0. We define a CFT with a boundary in EI such that the boundary condition is 
changed by a boundary operator. The random curve in H defined by SLE is growing 
so that it has states of one type to the left and of the other type to the right (the 
simplest way to view this is the lattice Ising model with the states defined as spin 
positions up or down). The mapping g satisfying ( 1271) 'unzips' the boundary. The 
primary operator that induces the boundary change with the conformal weight h 
is drift-less, and therefore, its expectation value does not change in time under the 
boundary unzipping. Hence all correlators computing with this operator remain 
invariant. Analogous considerations one may provide for the 'radial' version of SLE 
in the unit disk, slightly modifying the above statements. 
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